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Sympathetic cooling ol trapped ions has become an indispensable tool for quantum information 
processing and precision spectroscopy. In the simplest situation a single Doppler-cooled ion sym- 
pathetically cools another ion which typically has a different mass. We analytically investigate the 
effect of the mass ratio of such an ion crystal on the achievable temperature limit in the presence 
of external heating. As an example, we show that cooling of a single A1+ with Be"*", Mg"'' and 
Ca^ ions provides similar results for heating rates typically observed in ion traps, whereas cooling 
ions with a larger mass perform worse. Furthermore, we present numerical simulation results of 
the rethermalisation dynamics after a background gas collision for the Al^/Ca^ crystal for different 
cooling laser configurations. 



I. INTRODUCTION 



Sympathetic laser cooling of trapped ions is an impor- 
tant experimental tool in diverse fields, such as quantum 
information processing, chemistry and precision spec- 
troscopy. In quantum information processing, a sym- 
pathetic cooling ion can be used to cool the qubit ions 
in a multiplexed trap architecture [1], without affecting 
the quantum information store in the internal states of 
the qubit ions [H ^ . Sympathetic cooling of atomic and 
molecular ions enables studies of chemical reactions at 
cold and ultra-cold temperatures f4H6]. A major inter- 
est in sympathetic cooling of clouds of ions stems from 
spectroscopy of atomic [THH] and molecular ion species 
[TOl - fTS] with a complex internal level structure that can 
not be laser cooled directly. Precision spectroscopy of 
single or few spectroscopy ions sympathetically cooled 
by a well-controllable cooling or logic ion in a linear ion 
trap can be performed using quantum logic spectroscopy 
(QLS) [T11[T5]. Here, the laser-cooled logic ion not only 
provides sympathetical cooling, but also assists in the 
readout process of the spectroscopy ion after interro- 
gation. After its successful implementation in optical 
frequency standards based on aluminium ions |161 117j . 
several experiments including spectroscopy of molecular 
ions dHl IS], highly-charged ions [2DH12, super-heavy 
ions [23], and metal ions [21] using this technique have 
been proposed. One of the requirements for these appli- 
cations is efficient Doppler laser cooling of small linear 
ion crystals, which in its simplest form consists of two 
ions of (in general) unequal mass. Doppler laser cooling 
has first been experimentally demonstrated with trapped 
magnesium ions [5S] . It was soon realized that the strong 
mutual electro-static interaction between a laser cooled 
ion species and another species not interacting with the 
cooling laser, allowed sympathetic cooling of the latter 
in large clouds of ions [5^1 [17] . Crystallized linear chains 
of up to 15 ions have been sympathetically cooled with 
a single cooling ion of similar mass [28H30] . Sympathetic 
Doppler and even ground-state cooling of two-ion crys- 



tals with a mass ratio of up to three has been successfully 
implemented [TSl [3T] . Cooling of even larger mass ratios 
has been proposed by trapping the two species in sepa- 
rate potential wells [3ll|33]. The structure, dynamics and 
cooling of linear ion crystals composed of ions with un- 
equal mass have been investigated theoretically mostly 
in the context of quantum information processing |34l - 
I38j . with particular emphasis on the mode structure and 
the cooling rates. For applications, such as optical fre- 
quency standards requiring high spectroscopic accuracy, 
the lowest achievable temperature during Doppler cool- 
ing will determine relativistic and trapping-field induced 
shifts [T71 [55] . Motional heating of the ions in the trap 
due to electric field fluctuations [501 UHl HI] significantly 
modifies the achievable Doppler cooling temperature. 

Here, we develop an analytical model to study the 
achievable motional energy using sympathetic cooling in 
linear two-species two-ion crystals within the pseudopo- 
tential approximation in the presence of external mo- 
tional heating. The achievable cooling limit strongly de- 
pends on the strength of the external electric field fluctu- 
ations, the mass ratio between cooling and spectroscopy 
ions, and the trap parameters. We show that in particu- 
lar, the Doppler cooling temperature in radial direction 
is much more sensitive to a mass mismatch compared to 
the axial direction. We use this model to investigate the 
second-order Doppler shift for Doppler cooled ion crys- 
tals in an ^^Al+ quantum logic clock with different logic 
ion species. We show that next to the obvious choice 
of ^^Mg^-ions which excel because of their almost per- 
fect mass match, ^°Ca"'"-ions will perform similarly and 
in some situations even better than the Mg ions, owing 
to the lower achievable Doppler cooling temperature. Be- 
sides the cooling limit, the required cooling time after e.g. 
a collision with the hot background gas is an important 
aspect. We numerically simulate the cooling dynamics of 
a Ca+/Al+ ion crystal after such a collision event. Before 
crystallization of the two-ion crystal, the cooling rate is 
comparable to the single ion case. After crystallization, 
the motion of the ions is described in normal modes and 
the cooling rate is limited by modes that are only weakly 
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FIG. 1. Schematic of a hnear ion trap. The trap consists of 
four blade-shaped electrodes of which two opposing ones are 
connected to an RF-voltage while the other two are connected 
to ground. It also includes two endcap electrodes that are 
connected to a positive DC voltage. The line through the two 
endcaps defines the trap axis z and the other two axes (x, y) 
are chosen such that two blades lie on each axis. 



cooled by the logic ion. 

In section [IT] we briefly recall the low temperature dy- 
namics of the two-ion crystal and introduce analytic for- 
mulae for the six mode frequencies and the modal am- 
plitudes. In section [Tn] we expand the standard Doppler 
cooling model to the two-ion crystal case and investigate 
the effect of external heating on the temperature limit. In 



section IV we describe our cooling dynamics simulation 
and discuss the results before summarizing in section [V] 



II. NORMAL MODES OF A TWO-ION 
CRYSTAL 



approximation, a single ion of mass m and charge e in 
this trap will experience a potential of the form |45j 



TTt \ ^ 2 2,^ 22,^ 22 /r,\ 

U{x,y,z} = -muj^x + -mujyy + -mu^^z . (2) 



In this potential, the uJx,y,z denote the single ion trap fre- 
quencies, i.e. the frequencies of the ion's secular motion 
along the different axes. They are given by 



2eKUo 



id? 



y^w2 - (1 - a)ujl 



(3) 
(4) 
(5) 



where lo^, uJx and LOy are the axial and two radial trap 
frequencies, respectively, and 

describes the contribution of the RF potential to the ra- 
dial trap frequencies |43j . Introducing a factor e — iOp/ujz 
as in [3S| simplifies the radial trap frequencies to 



a to. 



(6) 



From the known trap frequencies, the a and e param- 
eters for a given ion with mass m = mi in a given trap 
can be derived. This allows the computation of frequen- 
cies for ions with different masses m2 in the same trap. 
While the axial trap frequencies simply scale with the 
square root of the mass, the radial trap frequencies ad- 
ditionally depend on e and a: 



In linear Paul traps the confinement of charged parti- 
cles is realised by two distinct electric fields [H] . The first 
is a rapidly oscillating 2D electric quadrupole in the ra- 
dial {xy) plane and the second is a static 3D quadrupole 
field providing confinement along the trap (z) axis. While 
the first is assumed to be purely radial and have no field 
components in the axial direction, the latter must have 
components in the radial directions additionally to its ax- 
ial field in order to satisfy Laplace's equation. Following 
the notation of [13] the total trap potential is given by 

y, 2, = y cos VlTt {l + 

+ .^o( ^''""^/'"^^' ) , (1) 

where Vq and Uq are the applied RF and DC voltages, 
R' K, R and d are the radial and axial trap dimensions, 
respectively, k is a geometric factor (44), < a < 1 
is a parameter indicating the radial asymmetry of the 
static field, and ^Ij- is the angular frequency of the rapidly 
oscillating RF-field (see Fig. [T]) . In the pseudopotential 



i^z.2 



l^x,2 




(7) 

(8) 

(9) 
1/2 to 



For the remainder of this work we assume a 
simplify the algebra. 

If two ions are simultaneously trapped in the same lin- 
ear Paul trap and strongly cooled to near K temper- 
ature, they will eventually crystallize at equilibrium po- 
sitions along the trap axis, equally spaced at a distance 

^0 = ( 32C/o7r£o ) ^^^P centre The 

remaining motion of ions 1 and 2 can then be described 
as small, coupled oscillations qi , q2 around these equilib- 
rium positions. Along every principal axis the motion 
consists of a superposition of an out-of-phase mode (o) 
where the two ions always move in opposite directions 
and an in-phase mode (i) where the two ions move in 
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the same direction. Following the approach of [3 
oscillations along a chosen direction are given by 



the 



qi (t) = Zibi sm{ujit + 00 + ^062 cos(woi + (t>o) (10) 
q2it) = sin{ujit + (pi) - cos{uiot + (fio) (11) 

where uJi^o,4'i,o are the angular eigenfrequencies and 
phases of the in-phase and out-of-phase modes, respec- 
tively, and bi^2 are the components of the normalised 
eigenvector of the in-phase mode, satisfying h\ + h\ — 1, 
in a coordinate system where the motion of the second ion 
is scaled by a factor of l/y/I with = 7712/7711. The Zi 



are the modal amplitudes (see Eq. (19), (20)). The calcu- 



lation of the modal frequencies and the 61^2-parameters 
can be performed similarily to for example [36] : For every 
ion the sum of the trap pseudopotential and the Coulomb 
potential due to repulsion from the other ion is developed 
around the equilibrium positions and the coupled equa- 
tions of motions are solved in lowest order, neglecting 
higher order non-linear couplings [201 135] . The results of 
this calculation for ions with different mass ratios ji and 
different e-parameters are given by: 



1 + ^ - ^Jl — ^ + fi'^ 



1 + A* + v^l - + 



ji + l_L^ — {1 + y?) — a 
/ + ^i"^ - e"^ [l + ^i'^) + a 

V V ' 

^i- + (-1 + +a 



2a 



(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



where the parameter 



\^,^-lf-2e^{^l-lY^l{l+^l)+^l^{l+{pi-l)^l) (i8) 



was introduced. The 62 parameters calculate as 62 = 
yjl — h\. Figure [2 shows the calculated eigenmode am- 
plitudes and frequencies in axial and radial directions. 
The heavier of the two ions has the largest amplitude 
for the mode with the lowest frequency, which is the ax- 
ial in-phase mode and radial out-of-phase mode. It is 
worthwhile noting that the radial mode amplitudes are 
much more sensitive to a change in the mass ratio. As 
a consequence, the radial motion of the ions is nearly 
decoupled for ion species with mass ratios /i < 0.25 or 
/i > 4 for typical traps, i.e. for one mode ion 1 has a 
large normal mode amplitude and ion 2 has a small one. 



whereas for the other mode the situation is reversed. The 
radial curves end at mass ratios where the radial out of 
phase mode frequencies reaches zero because for higher 
mass ratios the ion crystal turns from a linear axial con- 
figuration to a linear radial configuration [371 EH] ■ 

In this notation, the total energies of the two modes 
(i.e. kinetic -I- potential energy) along one direction 
(valid for radial and axial directions) are given by 



and 



(19) 
(20) 



III. COOLING LIMITS 

A. Doppler cooling of an isolated system 

In the following, we will briefly outline the usual deriva- 
tion of the Doppler cooling limit jl5M51] . before expand- 
ing the model to take into account an additional heat- 
ing rate in the next section. The model differs from the 
standard derivation in so far as it takes into account the 
modal structure of the two-ion crystal right from the be- 
ginning. 

Doppler cooling of all modes is achieved by exposing 
ion 1 to laser radiation of intensity / and an angular fre- 
quency uj = loq + A detuned by A from ion I's cooling 
transition with angular frequency wq, FWHM linewidth 
r and saturation intensity Iq. The fe- vector of the radi- 
ation is given by fc = {lx,ly,h)k = {lx,ly,lz)'-^o/c with 
the unit vector {lx,ly,lz)- 

Every scattering event of the radiation with the cooling 
ion will on average change the momentum of the system 
by hk. Due to the red detuning (A < 0) of the cooling 
light this is more likely to happen when the cooling ion 
moves towards the laser such that on average energy is re- 
moved from the ion crystal. The average energy loss rate 
can be calculated by averaging the momentum change 
over the periods of both oscillatory modes in each direc- 
tion. For sufficiently cold ions (|Za;fcgi| <C F) in the weak 
binding regime (F 3> Wi,o), these can (exemplarily in the 
x-direction) be shown to be (see Appendix IaI) 



dt cool 
dt cool 



2A/F 



/o(l + ///o + (2A/F)2)2 
/ 2A/F 
io(l + ///o + (2A/F)2): 



zjhluj (21) 



^^IbWo- (22) 



The competing heating rate (considering the statistical 
distribution of the momentum changes due to both the 
absorption and the spontaneous emission of the cooling 
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FIG. 2. Normal mode frequencies and normalized amplitudes for a two-ion two-species crystal. The square of the normalized 
amplitude (6i) and the in-phase [uJi) and out-of-phase (uJo) trap frequencies normalized to the axial trap frequency of a single 
ion of mass mi are shown for the axial and radial direction for different mass ratios ^ = 7712 /mi and e-parameters. 



photons) is given by (see Appendix [A|) 



^hcat ^ " To l2nH 1 + I/Io + (2A/r)2 



(23) 



dt heat ^ ^'^ - + ^ /o i2mi 1 + I/h + (2A/r)2 ■ 



(24) 



The cooling rate is a function of the ion's scattering 
rate and proportional to the square of the k-vector com- 
ponent (/2) of the cooling light along the considered mode 
direction, times the square of the motional amplitude of 
the cooling ion for that mode {zfb\, z^b"^). In contrast, 
the heating rate has a component from the isotropic pho- 
ton emission in addition to directed absorption from the 
cooling laser beam. 

The steady state solution is obtained from a balance 
between cooling and heating rates 



dE, 



dt heat dt cool 



(25) 



The resulting cooling limit can be expressed as 

?i(4A2 + r2(i + ///o))(i + 3;2) 



E, 



limit 



A8\A\n 



(26) 



for both modes. The cooling limit is independent of the 
mode eigenvectors 61^2 since the heating and cooling pro- 
cesses act in the same way on the motion of the ion crys- 
tal. As a consequence, modes with a small eigenvector 
component experience cooling and heating rates that are 
reduced by the same amount and will limit the time it 
takes to reach the steady-state temperature (see section 
IV). In the case of very low cooling intensity {I/Iq — )■ 0), 
optimum detuning (A — — r/2) and a laser that cools 
all directions symmetrically {l^ — ly — Iz = l/'v/S), this 
results in the well known Doppler cooling limit 



Ed = hr/2. 



(27) 



Very low cooling intensity leads to the lowest theo- 
retical cooling limit at the expense of long cooling times, 
since the cooling rate gets very small. In practice, cooling 
intensities close to the saturation intensity are typically 
used. A larger cooling rate renders the system more ro- 
bust if exposed to additional external heating. The max- 
imum cooling rate is achieved at / = 2/o and A = —r/2. 
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FIG. 3. Dependence of the coohng hmits on the angle of the 
cooling laser, fcy/fctotai denotes the component of the laser 
radiation parallel to the investigated axis. The other two 
axes are assumed to be cooled equally. 



However, this choice of parameters increases the cooling 
limit by a factor 2. 

Aligning the cooling beam along the direction of a par- 
ticular set of modes, cooling below the Doppler-cooling 
limit given by Eq. ( [27| in this direction is possible. How- 
ever, the cooling limit in the other directions will be 
strongly increased as can be seen in Fig. |3] and has been 
discussed in [52]. This is a direct consequence of the 
cooling rate of a specific mode being dependent only on 
the k-vector projection along this direction, whereas the 
heating rate has a contribution from the isotropic spon- 
taneous emission heating. 

Precision spectroscopy and atomic frequency stan- 
dards based on trapped ions require small kinetic ener- 
gies to achieve small second-order Doppler shifts, which is 
particularly important for light ions such as ^^Al+. This 
shift is given by [S3] 



V 



2c2 



(28) 



where (w^) is the average value of the square of the ion's 
velocity and c is the speed of light. This relation allows a 
direct mapping between cooling limits and corresponding 
relativistic Doppler shifts. In the absence of external 
heating rates, both ions will have the same energy given 
by Eq. (27) and the second-order Doppler shift of the 



clock ion is 



At/ 

V 



7712 C 



(29) 



In the case of ^^Al+ sympathetically cooled by ^°Ca+, 
this results in shifts of 9.2 x 10~^^ for each of the three 
directions. For the radial modes, this number has to be 
multiplied by a factor of roughly 2 (see Appendix |b| to 
account for the kinetic energy in the intrinsic micromo- 
tion of the radial motion of the clock ion which follows 
from generalisation of the results of [33] . The only influ- 
ence the selected cooling ion species has on this result is 



the linewidth of its cooling transition. Furthermore, the 
mass ratio determines the time it takes to reach steady 
state, but has no influence on the cooling limit. This 
changes as soon as external heating rates are included in 
the model. 



B. Doppler cooling with external heating rates 

External heating rates of cooled ion crystals are as- 
sumed to be mostly due to stochastic electric field fluc- 
tuations [HJ 1331 [S5l [SB]. If the characteristic distance 
between the electrodes and the ions is much larger than 
the distance between the ions in the crystal, the electric 
field across the ion crystal can be assumed constant and 
the heating rate due to field fluctuations can be written 
as .36' : 



dEi 

dt fluct 


q'SE 
Ami 




1 


dEo 

dt fluct 


Ami 




1 



bi 



(30) 
(31) 



Here Se = Se{(jj) denotes the electric field spectral 
density, which is assumed to be spectrally constant for 
the relevant Wi_o- The homogeneous field fiuctuations 
will only couple to centre of mass motion and therefore 
much more strongly to the in-phase mode than the out- 
of-phase mode. In fact for /it = 1 it follows that bi = 62 
and the out-of-phase mode is not heated at all. The to- 
tal energy injected in one direction of motion is obtained 
by adding the in-phase and out-of-phase heating rates, 
which turns out to be proportional to m^^ -|-m^^, show- 
ing the advantage of heavy cooling and clock ions. 

The cooling limit in the presence of external heating 
is obtained by including the heating rates Eq. (30) and 
Eq. (31) in the steady state condition Eq. ([25[). The 
steady state energy limit for the in-phase mode (substi- 
tute 61 by 62 and 62 by —61 for the out-of-phase-mode) 
is given by: 



E, 



limit, 1 



r (1 + (2A/r)2 + i/Iq 

48|A|W//oZ2fc2 



rh^i/io{3il + i)P 



1 



2(^1 



1 



b2fq^SE 



(32) 



(2A/r)2 + I/Io 
The additional heating modifies the steady state solution 



for Doppler cooling (Eq. (26)) by breaking the symme- 



try between cooling and heating processes: The photon- 
induced heating and cooling rates have the same depen- 
dence on the mode amplitudes, resulting in a Doppler 
limit independent of this parameter. The external heat- 
ing rate has no cooling component and a more complex 
dependence on the mode amplitudes, resulting in a cool- 
ing limit that is a sensitive function of the modal ampli- 
tudes and therefore the e and fi parameters. As a result. 
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modes where the cooling ion has a large relative ampli- 
tude are cooled more efficiently than modes where it has 
a small amplitude. 



The two terms in the square bracket in Eq. ( 32 ) give 
the contribution of the photon and the external heating 
rate to the total heating rate, respectively. It is instruc- 
tive to investigate the two extreme cases, in which either 
one of the contributions dominate. Defining an electric 
field spectral density Seo = Th^k'^/I2q^ [57] for which 
the two heating rates in a symmetric {jj. = 1) ion crystal 
and symmetric cooling in all three directions {Ix — 
equal ^,„,, = ^,,,,(A = -r/2,/ ^ 21,), the two 



dE, 

dt fiuct 

cases are 



1. Se ^ Seo- In this case the photon heating rate at 
the point of maximum cooling is much larger than 
the external heating and the latter can therefore be 
neglected. By choosing the intensity < ///g < 2 
optimally, the cooling limit of Eq. ( 26 ) is recovered: 



;ir/2 < Eunnt < hT 



(33) 



Se ^ Seo- In this regime the photon heating rate 
can be neglected and the cooling limit will pro- 
portionally depend upon the electric field spectral 
density and the mode amplitudes. Here the opti- 
mum intensity is that of maximum cooling rates 
(/ = 2/o). 



^limit,^ 



-^limit,o 



q'SEih + ^b2fT (1 + (2A/r)2 + I/IoY 
16bl\A\hllk^I/Io 

q'SE{b2 - ^b.rr (1 + (2A/r)2 + I/I, 



(34) 

2 



l6bl\A\nP^k^I/Io 



(35) 



In both regimes the best cooling performance is always 
achieved at a detuning of A = —r/2. In the remainder of 
this paper we will look at external heating rates of up to 
Se = 0.2Seo that are typical for macroscopic ion traps. 

For the evaluation of the second-order Doppler shift 
due to motion along one spatial direction, it is not the 
total kinetic energy in a certain mode that is of interest, 
but rather the total kinetic energy in the secular motion 
of the clock/spectroscopy ion. This energy is given by 
the sum of the energies in both modes along that spatial 
axis, weighed by the relative fraction of clock ion energy 
to the total energy in the respective mode: 



{blE, 



b{Eo 



(36) 



Figure |4] shows a plot of E^ as a function of the mass 
ratio ji for axial motion. The y-axis is normalised to 
the energy of the clock ion at the Doppler-cooling limit 
without external heating and the electric field spectral 
densities Se are given as multiples of Seo- In this fig- 
ure, the intensity was optimised for each data point to 
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FIG. 4. Normalized axial clock ion energy plotted against the 
mass ratio ^ of the ion crystal. The energy Ec of the clock 
ion is the sum of the clock ion energy in both modes. The 
calculations were performed for varying electric field spectral 
densities Se- The cooling laser intensity I / lo has been opti- 
mized for each value of Se- The dashed lines show different 
logic ion species for an Al""" clock. 



achieve the lowest energy in the clock ion. The minimal 
clock ion energy for different electric field heating rates is 
always achieved at a mass ratio of jj. = 8/11. This plot is 
independent of F and the actual values of mi, m2- How- 
ever, the normalization factor and therefore the absolute 
value for the cooling limit depends on the linewidth of 
the cooling ion. This doubles all cooling limits for Mg"*" 
with respect to all other ions, owing to its twice as large 
cooling transition linewidth compared to the other ions. 
The figure shows that the axial clock ion kinetic energy 
is only a weak function of the mass ratio. This is a direct 
consequence of the weak dependence of the axial mode 
amplitudes on this parameter, as shown in Fig. [2j 

Figure [5] shows the normalized clock ion energies in one 
of the radial directions for varying mass ratios and dif- 
ferent e parameters, assuming a radially symmetric trap 
[a = 1/2). The stated e is always that of a single clock 
ion in the trap. The shaded regions in the graphs show 
the areas in which the crystal is not stable (see HI I and 
hence can be ignored. The graphs show that mass ratios 
slightly larger than 1 always result in the lowest possible 



clock ion kinetic energy. For a given mass ratio, the en- 
ergy of the clock ion in the presence of external heating is 
lowest when operating the trap close to instability of the 
linear axial configuration. The reason for this is that the 
larger the radial confinement, the more the modal am- 
plitude in the radial modes tend towards &i w 1, 62 ~ 
or vice versa (see Eq. (17) and Fig. [2|. In that case, the 
mode with negligible motion of the cooling ion will cool 
very poorly, leading to an elevated steady state tempera- 
ture in the presence of external fields. In this regime, the 
comparison of Fig. |4] and [5] (the y-axes are normalised to 
the same energies) shows that the heating-induced radial 
clock ion energy is larger than the axial clock ion energy 
and will therefore dominate the clock frequency shifts. 
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FIG. 5. Normalized radial clock ion energy plotted against the mass ratio fi of the ion crystal for radial modes. The energy Ec 
of the clock ion is the sum of the clock ion energy in both radial modes. The shaded region indicates parameters for which the 
ion crystal is no longer linear in axial direction. The cooling laser intensity I/Io been optimized for each value of Se- 

Choosing an e parameter close to the instability regime 
improves the cooling limits, but at the same time in- 
creases the intrinsic micromotion amplitude of the radial 
modes (see Appendix [b|) . Fig. [6] shows the normalized 
total kinetic energy (equal to the sum of secular and 
micromotion kinetic energy) of the clock ion in an ion 
crystal that has been cooled to the Doppler cooling limit 
(i.e. without external heating) plotted against the crystal 
mass ratio. 

Since the micromotion kinetic energy contributes to 
the second-order Doppler shift in the same way the sec- 
ular energy does, it is not advisable to perform spec- 
troscopy in a trap operating close to the unstable regime. 
The operation in a trap with very strong radial con- 
finement is equally bad because of the inefficient cool- 
ing of the weakly damped radial mode. For best cooling 
performance it is therefore advisable to optimise the e- 
parameter according to the observed heating rate in the 
trap. 

A comparison of the total second-order Doppler shift 
(including the scaling factors due to micromotion for the 
radial modes) in a linear ion trap for Al+/X+-ion crystals 
with different possible cooling ions X"*" is shown in Fig. 
[7j The y-axis gives the normalised electric field spectral 
density SE/SEo,Ca [53 . A value of Se/Seo = 0.02 cor- 
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FIG. 6. Total clock ion energy (secular energy -I- micromotion 
energy) of a clock ion in a Doppler cooled crystal along one 
radial direction. Here the mass ratio ^ and the e-parameters 
were varied and the energy was normalised to the Doppler 
cooling energy Eu- The absence of external heating was as- 
sumed. Due to large micromotion contributions the total en- 
ergy diverges at the points where the crystal becomes instable. 



responds to a radial heating rate of a single Ca+ ion in a 
trap with radial trap frequencies of 2.5 MHz of roughly 
1500 quanta per second and can be regarded as an upper 
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limit for most ion traps used for spectroscopy. The 
value as well as the Doppler laser detuning A and inten- 
sity I/Iq were optimised for each value of Se- 

The graph shows that an Al+/Mg+ crystal suffers the 
least from external heating since its mass ratio is very 
close to 1. However, for traps with fairly low heat- 
ing rates, the Doppler shift in Al+/Ca+ and Al+ZBe'^ 
traps will be lower than that of Mg^-systems because 
the linewidth of the cooling transition of Mg^-ions is ap- 
proximately a factor of 2 larger than that of Ca+ and 
Be+-ions. The two other species Yb"*" and Sr+ suffer 
from their small mass ratios and require traps with very 
small heating rates to reach comparable performance. 

At this point it should be noted, that not the abso- 
lute value of the Doppler shift, but rather its uncertainty 
is relevant for the clock performance. However, since 
this uncertainty can be assumed to scale with the abso- 
lute shift, a reduction of the latter will result in a re- 
duction of the former and therefore lead to better clock 
performance. If Doppler cooling is turned off during spec- 
troscopy, careful modelling or measurement of the heat- 
ing process is required to establish the uncertainty. 



IV. COOLING TIMES 

An atomic clock will reach it's maximum performance 
in terms of stability, if the clock transition is probed with- 
out any dead time between consecutive readings 58 60j . 
One contribution to the dead time is the time it takes to 
cool the ion crystal in case a collision with a hot back- 
ground gas particle occurred. The large energy transfer 
to the ion crystal leads to typical temperatures on the 
order of room temperature, resulting in the decrystal- 
lization of the ion crystal. As a consequence, it is neces- 
sary to take the non-linear contributions of the Coulomb 
forces into account to evaluate the cooling dynamics at 
these high temperatures. The nonlinearities are usually 
neglected in the small oscillation approximation used to 
describe the low temperature dynamics of the system 
that we have used in the previous sections. The dynamics 
of small ion crystals during laser cooling in various tem- 
perature regimes have first been studied in the context 
of phase-transitions and chaos theory [HI [62] ■ 

Here, we numerically solved the equations of motion 
resulting from the complete ponderomotive two ion po- 
tential 

V{xi,yi,zi,X2,y2,Z2) = ^TOi (w^^xi + cj^i2/i +^lizl) 



47reo ^{xi - X2Y + {yi - y2)2 + (zi - Z2)^ ' 



(37) 



expect only minor modifications to the results presented 
here when micromotion is included in the treatment j63j . 
We used an adaptive Runge-Kutta method where the 
step size was reduced significantly every time the ions 
got close enough to explore the high non-linearities of the 
1/r-Coulomb potential. Cooling was incorporated into 
the calculation by multiplying the instantaneous scatter- 
ing rate R at every time step with the instantaneous step 
size dt at the same step and comparing the resulting value 
with a random number rn(0, 1) between and 1 from a 
"Mersenne Twister" -type pseudo random number gener- 
ator [53] . A scattering event that changed the velocity 
Vi of the cooling ion according to the momentum change 
due to both absorption and spontaneous emission was 
therefore calculated whenever the condition 



Rdt = 



IJIo 

2l + ///o + (2(A-fc.t;i)/r) 



^dt < rn(0, 1). 



where cukj is the trap frequency of ion j in direction k G 
x,y,z. and kj are the position coordinates. For the sim- 
ulations, the RF potential was neglected. However, we 



(38) 
was fulfilled. 

With a small enough step size the simulation turned 
out to be very robust and reproduced two-ion cooling lim- 
its as well as theoretical curves for the cooling duration 
of single ions (e.g. [BS]). 

A typical collision event with the background gas will 
be an elastic collision of a hydrogen molecule with either 
one of the two cold ions in the crystal. The maximum en- 
ergy that can be transmitted in an elastic collision of two 
particles with masses mi, m2 if one of them is initially at 
rest is given by times the kinetic energy of the 

hot particle. In case of the collision of a cold aluminium 
ion and a hydrogen molecule at 300 K this means that 
a maximum energy of « 0.26 x 3/2kB x 300 K can be 
transferred to the aluminium ion. Here ks denotes the 
Boltzmann constant. This corresponds to a temperature 
of the ion crystal of « 19.3 K. Cooling an ion crystal 
from these high temperatures is aided by adding a far 
detuned laser beam additionally to the standard cooling 
beam at A = — r/2. This enhances the scattering rate of 
the rapidly moving cooling ion with large Doppler detun- 
ing. The optimum detuning of the second laser was found 
by simulating the cooling time for an aluminium/calcium 
ion pair starting at a temperature of 19.3 K as a function 
of the detuning. The results are shown in Fig. [Sj 

Both lasers were assumed to be directed onto the ion 
from the (1, 1, 1) direction and both had one saturation 
intensity at the position of the ion. To compare the per- 
formance of the different detunings, the time to reach a 
crystal energy of 1000 Doppler cooling limits was simu- 
lated. This energy roughly corresponds to the energy at 
which the ions crystallize. Fastest cooling was achieved 
at a detuning of ^-100 MHz. This optimum detuning of 
the second laser beam proved to be very robust against 
changes in the cooling parameters, such as a change in 
linewidth of the cooling transition, the intensity of the 
cooling laser, or the mass of the cooling ion. Since 
these changes mostly affect the scattering rate, the times 
needed to cool vary strongly. The general characteris- 
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FIG. 7. Comparison of secular motion-induced second-order Doppler shifts for Al^/X^-crystals with different cooling ions 
={Be^, Mg"^, Ca^, Sr+, Yb"*"} in the presence of external heating with an electric field spectral density of Se- The right 
Figure shows the parameters e and I/Iq for best cooling performance for an Al"'"/Ca'*'-crystal. Those optimum parmaters for 
a spectral energy density of S = O.OOSSbo are given by 2.63,0.40,3.17,0.05,1.12,0.87 and 1.04, 1.61 for the Be, Mg, Sr and Yb 
crystals, respectively. 




50 100 150 200 250 300 350 400 450 500 

Detuning of the second laser [MHz] 

FIG. 8. Time needed to cool an Al^/Ca'*"-ion pair from 19.3 K to 1000 Doppler limits. The y-axis shows the detuning of a 
second laser at saturation intensity that is cooling the ion crystal in addition to the standard — r/2-detuned beam. 



tics of the curve with the minimum close to -100 MHz, 
however, stays the same. This can be explained by the 
fact that in a well-thermalized regime, the cooling ion 
will carry half the crystal energy and the optimal de- 
tuning should only depend on the average cooling ion 
velocity in the direction of the incoming laser. This is in- 
dependent of the transition lincwidth or the intensity of 
the cooling light. The change of mass of the cooling ion 
to smaller values shifts the optimal detuning to slightly 



higher values. This is because a lighter cooling ion has 
higher average velocities if it has the same energy as a 
heavier ion. 

The difference in cooling dynamics between a two- 
and a single-ion system is shown in Fig. [9j Besides the 
Al+ZCa"*" pair, the results for a single Ca+ ion in the trap 
with the second laser detuned by 140 MHz is shown. This 
slightly higher optimum detuning arises, since a single 
ion will on average have higher velocities than the cool- 
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ing ion in a two-ion system since the latter can exchange 
energy with the clock ion. For the simulation, the start- 
ing energy of the single Ca+ ion was set to the same value 
as that of the two-ion pair, although a collision with a 
hydrogen molecule would lead to a lower initial temper- 
ature. 

The cooling rates match extraordinarily well at the 
beginning of the cooling process. This is because in the 
high temperature regime the two ions in the crystal col- 
lide very often so that all motional modes are thermal- 
ized almost instantly. Therefore no poorly cooled modes 
exist that would decrease the cooling rates. However, 
as soon as the ion pair crystallises (around 7 ms), no 
collisions thermalize the modes and the cooling rate de- 
creases significantly owing to the weakly damped radial 
modes (see inset in Figure [9| . The crystallisation effect 
can be visualised by plotting the energy of the Ca"*" and 
Al"*" separately over the cooling time, as shown in the 
inset of Figure [9j Once crystallisation occurs, the en- 
ergy exchange between the motional modes of the two 
ions is prohibited and the cooling rate is determined by 
the mode amplitudes of the cooling ion (for a theoreti- 
cal description of the crystallisation see for example [52] 
for a measurement of the damping of the weakly damped 
modes see the supplementary material of |16j). 

Summarizing the simulation results of Fig.|8]and[9j an 
optimum cooling strategy after a background gas colli- 
sion can be derived: The fastest cooling is achieved by 
ramping the laser detuning according to the actual en- 
ergy of the ion crystal, such that maximal scattering rates 
are maintained. Additionally, the intensity of the cooling 
laser should be large to increase the rates even further. 
Once crystallization occurs, the laser should be ramped 
to a detuning and an intensity for which the lowest ener- 
gies are achieved (typically A«r/2, / ^/q). 

In our simulations the time needed to cool an 
Al+/Ca+-crystal to 2 Dopplcr-Cooling limits is roughly 
90 ms (which corresponds well to previous studies [66' 
while the time it takes to reach crystallisation is only 
7 ms. A possibility to reduce the final cooling time is to 
couple the well-damped to the poorly-damped modes by 
a static radial electric field as has been demonstrated in 
|16j . Even without such a scheme, a Doppler cooling time 
of < 100 ms corresponds to less then one measurement 
cycle in state of the art quantum logic clocks (i.e. cy- 
cle time 230 ms TT) and since such a random collision 
should only occur roughly every 100 s in a good vacuum 
environment, its contribution to the total dead time will 
be insignificant compared to dead times due to detection 
and state preparation. However, it seems necessary that 
the collision events are detected and actions are taken 
to ensure fast recrystallization, since the typical Doppler 
cooling time used in single ion experiments are on the or- 
der of ms, which will not be sufficient to cool the crystal 
to its cooling limit. 



V. SUMMARY 

In this paper we examined the steady state sympa- 
thetic cooling limits with and without external heating 
of different two-ion two-species ion crystals with an em- 
phasis on second-order Doppler shifts that pose a limit 
to precision spectroscopy. We find that in the absence of 
external heating the Doppler cooling limit can be reached 
in all six normal modes of the ion crystal, independent 
of the mass ratio between the two ions. However, with 
additional external heating, the cooling limit becomes a 
sensitive function of the normal mode amplitudes. The 
modes where the cooling ion has a large amplitude are 
cooled most efficiently and reach a low cooling tempera- 
ture. We find that operating the trap close to instability 
with respect to a flip of the linear-axial to the linear- 
radial orientation enhances the amplitudes of the critical 
radial modes. This effect has to be balanced with the in- 
creased micromotion in this regime. Specifically, we find 
that three logic ion species candidates, ^^Mg+, ^Be"^ and 
'*''Ca+, are most suitable as sympathetic cooling ions for 
an aluminium ion clock, even in the presence of moderate 
external heating rates. The former is a good choice be- 
cause its mass ratio relative to aluminium is close to one 
which means that no poorly-damped modes exist that 
would make the crystal vulnerable to high external heat- 
ing rates. The latter two have a cooling transition with 
a smaller linewidth than the ^^Mg+-transition resulting 
in lower Doppler-limits and therefore lower second-order 
Doppler shifts. For traps with very low heating rates, 
even heavier sympathetic cooling ions, such as ^®Sr+ will 
perform well. Furthermore, we investigated the cooling 
time of the "^^Ca+Z^^Al^ unequal mass ion crystal af- 
ter a collision with background gas. Monte-Carlo cooling 
simulations taking into account the ponderomotive trap 
potentials and the Coulomb-potentials revealed cooling 
times of roughly 100 ms, which would not pose a limita- 
tion for the clock stability. However, these events should 
be detected immediately and counteracted by applying 
a second (further detuned) laser beam for rapid recrys- 
tallization. We therefore believe that an ^°Ca+/^'^Al+ 
quantum logic clock will lead to a similar clock perfor- 
mance as 25Mg+/27Al+ clocks and might even outper- 
form them in terms of second-order Doppler shifts for ion- 
traps with low heating rates. Furthermore, the presented 
results are relevant for high precision spectroscopy of 
other sympathetically cooled ion species, such as highly- 
charged or molecular ions. 



Appendix A: Derivation of the Cooling and Heating 
Rates 



The cooling rate of a two-ion crystal mode can be cal- 
culated by determining the energy change in this mode 
for every scattering event. According to Eq. (19) the 



energy change in an in-phase mode (analogously for out- 
of-phase modes) with modal amplitudes Zi before and 
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Time after beginning of Cooling [s] 



FIG. 9. Cooling evolution of an Al^/Ca^ ion pair after a collision event with an H2-molecule (blue line). The red line shows 
the cooling evolution of a single Ca^ ion in the same trap with the same initial energy. The horizontal dotted lines denote 
the Doppler cooling limit and the inset shows the energies of the calcium and aluminium ion before and after crystallisation 
occured at around 7 ms. 



after the absorption of a cooling photon is given by 
AE, = (zf - zf) . 



(Al) 



zl can be calculated by expressing the modal amplitude 
Zi as a function of the positions qi , q2 and velocities vi = 
QiiV2 = 92 along the axis of the relevant modes of the 
two ions before the absorption 



^ ^ iVJj'b2V2 + biVi)'^ + Ujf{y/Jlb2Z2 + biziy ^^^^ 



and then adding the velocity change dv = Ir^hk/mi (ex- 
emplarily for the a;-direction) to the cooling ion velocity 

, _ / (V^^2t'2 + bi{vi + dv))^ + Ujf{^JIb2Z2 + biZiY 



(A3) 



Substituting these equations into Eq. (All AiJ^ gives 
1 



AEi = -TOi6]^dw +mi&idw(6iWi + \/Jib2V2) 
1 



(A4) 



-mifej^dw + rriibidvuJiZi cos (ujit + (pi). (A5) 



The first term of this equation is a constant heating 
that is taken care of in the heating rate. The sec- 
ond term corresponds to mode cooling if cos {lo it + (pi) 
is negative and heating if it is positive. The rate R 
at which scattering events occur is given by Eq. ( 38 1 . 



The cooling rate is computed by averaging the product 
R X {/S.Ei — }^mib\dv^) over the oscillation periods of all 



six crystal modes 
dE, 

dt Cooling 

277 

1 



{2n) 



R{AEi - lmibldv^)d^(f>i d^po- (A6) 



The integration over all six modes is necessary because 
the scattering rate depends on the product k ■ Vx and 
therefore on the cooling ion velocity in all three spatial 
dimensions. This integral cannot be solved in general but 
for cold crystals, where \kvi | ^ T in all three dimensions, 
the integrand can be expanded and higher order terms 
in kvi/r can be neglected, leading to Eq. (21 1. 

The heating rate computes similarly by averaging the 
product R X ^mibidv^ over the six modes which leads in 
first order to a velocity independent term 



r 



lo 12mi 1 + I/Io + (2A/r)2 



(A7) 



To this heating due to the absorption of the photons one 
has to add a contribution of the spontaneous emission. 
This is exactly as large as the absorption effect but does 
not depend on the direction of the cooling laser. Assum- 
ing an isotropic emission it is given by 



I 



bl 



lo 12mi 1 + I/Io + (2A/r)2 



(A8) 



so that the sum of Eq. (ATI and (AS) gives the total 
heating rate as in Eq. ( 23 ) . 
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Appendix B: Doppler Shift Contribution of Intrinsic 
Micromotion 



To calculate the amplitude of the micromotion in the 
radial modes of the two-ion crystal, we generalize the 
derivation for single ions given in |54j . The force on a 
single ion in radial direction (we exemplarily us the x- 
direction in the following) is given by 



eC/n 



x^^cosfitt (Bl) 



Fx = -d^e^ix, y, z, t) = 2xa ^ 

where ?/, 2;, i) from Eq. ([ij was used. In two-ion 
crystals an additional force arises due to the Coulomb 
potential Uq between the two ions 



Uc{xi,X2,vi,y2-, zi, Z2) = 

^ 1 

47reo ^{xi - 2:2)2 + {yi - yz)^ + (zi - Za)^ ' 

The Coulomb-force in x-direction is given by 

Fc,x = ~dxUc{xi,X2,yi,y2,zi,z2) = 
e2 1 



(B2) 



47re 



y/{xi - X2Y (2/1 - 2/2)2 _^ - Z2)2 



;\xx-X2). 



(B3) 



We now assume that the ion crystal is cold enough so 
that the ions oscillate with small amplitudes around their 

equilibrium positions xf^ = (0,0, zq), '^^^ = (0,0, —zg) 
/ 2 \ 1/3 

where zq — ( 3^^^^ j (see for example 03). Linearis- 



ing of Eq. (B3) in {x\ — X2) around aj^^^ajj"^ gives 

\ . . eUo 



F, 



C,x 



r(xi 



X2] 



(P 



If we now consider only a single mode (i.e. the in-phase 
or the out-of-phase mode), we can write the x-motion of 
the second ion as a constant factor /3 times the motion 
of the first ion 



X2 — (3xi 



(B5) 



where f3 = b2/{y/Jj.bi) for in-phase and (3 = —61/(^/162) 
for out-of-phase modes as follows from Ec^. ([lO| and ([TT|) . 



Replacing X2 in Eq. (B4) and adding the Coulomb- 



force to the force due to trap potentials, the total force 
Gx on the ion is given by 



Gx — Fx + F( 



c,x=x,^{2a + l-l3) 



Introducing the Ux 



Xij^ COSiltt. 

(B6) 

aeUp „„J „ _ 2eVo • 
—pTFrr anu qx — rniR^Sl^ 



[54j . the equation of motion for xi can be given in typical 
Mathieu form 



Xi 



1 



1-/3 
2a 



2qx cos (r^T^) ) ^^^1 ~ 0- 



(B7) 



Compared to the equation of motion for a single ion 

{xi + {qx + 2qx cos (Oyi)) -^Xi = 0), the only difference 
is that the Qx is varied by a factor (1-|-(1 — /3)/(2a)) in the 
case of a two-ion crystal mode. Accordingly the first or- 
der solution [55] (valid for < 1, a2,(l-|-(l-^)/(2a)) < 
1) of the single ion equation is still valid, if that factor is 
included: 



xi{t) « Xih smujx,it + (l)x,i (1 + y cos{D,Tt)j (B8) 



In this equation the amplitude Xibi of the oscillation was 
already chosen to comply with Eq. ( 10 1. ujx,i denotes the 
in-phase mode frequency and is given by 



axil + {I -P)/ {2a)) 



(B9) 



The average squared velocity of this motion is given by 



1 



2' 9^+2^.(1 + ^5^)^ 



1 



{x,bi)''uj%, ( 1 + 



2e2 



2e2 - 2a - (1 - /3) 



(BIO) 



where in the last step the qx and ax were substituted 



{X1-X2). (B4) by the = 



factor that is used throughout this 



paper. Equal ion [13 1 1 shows that for the radial modes the 
average squared velocity is not just given by the secular 
motion but that a micromotion term of the same order 
of magnitude has to be considered as well. For strong 
radial confinement ^ 1, the micromotion contribution 
approximately equals the contribution from the secular 
motion. However, for weak confinement, it can become 
much larger. 

Figure [T0| shows a plot of the relative kinetic energy of 
the micromotion of the clock ion compared to its secu- 
lar energy for one pair of radial modes and for different 
e-parameters (again the stated e is that of a single clock 
ion) in the absence of external heating. The graph shows 
that the out-of-phase fractional micromotion energy con- 
tribution can get very large, when the trap is operated 
close to the instability limit. 

The micromotion contribution factor for the second 
ion (which is in this paper the clock ion and therefore 
the interesting one), can be calculated identically but 
the £2 = \/ mi/TO2e and /32 — 1//3 parameters for the 
second ion have to be used. Expressed in terms of the 
e and 61,62 parameters, the average squared velocity of 
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FIG. 10. Relative micromotion energy of the clock ion in the in-phase and out-of-phase radial modes. The energy is normalised 
to the secular energy of the clock ion in the respective mode. In this graph the absence of external heating was assumed. The 
dashed lines denote the asymptotic behaviour for very large radial confinement (e — > oo). 

the clock ion is therefore given by 

•2 1 / ^ ^2 2 2eV^ \ 

< X,, > = -(..5.) -.{l+2,2/^_2Ml-^/AIfel/62)j 

(BU) 

< > = -i.M) -..(;+2,2/^_2Ml+VM^2/6t)j 

(B12) 

for the in-phase and out-of-phase mode, respectively. 
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